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We provide a systematic construction for all n-particle 
form factors of the SU(N) 2 /U(l) N - 1 -homo geneous Sine- 
Gordon model in terms of general determinant formulae for a 
huge class of local operators. The ultraviolet limit is carried 
out and the corresponding Virasoro central charge together 
with the conformal dimensions of various operators are iden- 
tified. The renormalization group flow is studied and we find 
a precise rule, depending on the relative order of magnitude 
of the resonance parameters, according to which the theory 
decouples into new cosets along the flow. 

PACS numbers: ll.lOHi, ll.lOKk, 11.30Er, 05.70.Jk 



I. INTRODUCTION 

For most integrable quantum field theories in 1+1 
space-time dimensions it remains an open challenge to 
complete the entire bootstrap program, i.e. to com- 
pute the exact on-shell S-matrix, closed formulae for the 
n-particle form factors, identify the entire local opera- 
tor content and in particular thereafter to compute the 
related correlation functions. Recently we investigated 
§-| a class of models, the S , t/(3) 2 /t/(l) 2 -homo geneous 
Sine-Gordon model Q] (HSG), for which this task was 
completed to a large extend. In particular we provided 
general formulae for the n-particle form factors related 
to a huge class of local operators. In order to understand 
the generic group theoretical structure of the n-particle 
form factor expressions it is highly desirable to extend 
that analysis to higher rank as well as to higher level. 
One of the main purposes of this manuscript is to do the 
former, that is the investigation of the SU(N) 2 /U(1) N ~ 1 
case. This model may be viewed as the perturbation of 
a gauged WZNW-coset with Virasoro central charge 
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by an operator of conformal dimension A = N/(N + 2). 
The theory possesses already a fairly rich particle con- 
tent, namely N — 1 asymptotically stable particles char- 
acterized by a mass scale rrii and N — 2 unstable par- 
ticles whose energy scale is characterized by the reso- 
nance parameters <xy (1 < i,j < N — 1). We relate the 
stable particles in a one-to-one fashion to the vertices 
of the SU ( A)-Dynkin diagram and associate to the link 
between vertex i and j the N — 2 linearly independent 
resonance parameters ov,- . 



We find that once an unstable particle becomes ex- 
tremely heavy the original coset decouples into a direct 
product of two cosets different from the original one 
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Equivalently we may summarize the flow along the 
renormalization group trajectory with increasing RG- 
parameter r to cutting the related Dynkin diagrams at 
decreasing values of the ct's. For instance taking c^+i to 
be the largest resonance parameter at some energy scale 
the following cut takes place: 
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Using the usual expressions for the coset central charge 
the decoupled system has the central charge 
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Our manuscript is organized as follows: In section II 
we present the main characteristics of the HSG-scattering 
matrix. In section III we systematically construct so- 
lutions to the form factor consistency equations, which 
correspond to a huge class of local operators. In section 
IV we investigate the renormalization group flow of the 
Virasoro central charge, reproducing the decoupling (|^). 
In section V we compute the operator content in terms of 
primary fields of the underlying conformal field theory. 
In section VI we investigate the RG flow of conformal 
dimensions. Our conclusions are stated in section VII. 



II. THE S-MATRIX 

The prerequisite for the computation of form factors 
and correlation functions thereafter is the knowledge of 
the exact scattering matrix. The two-particle S-matrix 
describing the scattering of two stable particles of type i 
and j, with I < i, j < N — 1, as a function of the rapidity 
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9 related to this model was proposed in H. Adopted to 
a slightly different notation it may be written as 



F 0\tm(l+,l-), 



S ij (6) = (-lY 



tanh - I 9 + cr. 
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The incidence matrix of the S l £/(-/V)-Dynkin diagram is 
denoted by /. The parity breaking which is characteris- 
tic for the HSG models and manifests itself by the fact 
that Sij ^ Sji, takes place through the resonance pa- 
rameters Gij = —<Jji and the colour value Cj. The latter 
quantity arises from a partition of the Dynkin diagram 
into two disjoint sets, which we refer to as "+" and "— " . 
We then associate the values Ci = ±1 to the vertices i 
of the Dynkin diagram of SU (N) , in such a way that no 
two vertices related to the same set are linked together. 
Likewise we could simply divide the particles into odd 
and even, however, such a division would be specific to 
SU(N) and the bi-colouration just outlined admits a gen- 
eralization to other groups as well. The resonance poles 
in Sij(6) at (9n)ij — —<Jij — in/2 are associated in the 
usual Breit-Wigner fashion to the N — 2 unstable parti- 
cles as explained for instance in |t|||[|. It is important 
for us to recall that the mass of the unstable particle 
Mg formed in the scattering between the stable particles 
i and j behaves as Mg ~ ek^l/ 2 . There are no poles 
present on the imaginary axis, which indicates that no 
stable bound states may be formed. 

It is clear from the expression of the scattering ma- 
trix (|]), that whenever a resonance parameter Cij with 
Iij 7^ goes to infinity, we may view the whole system 
as consisting out of two sets of particles which only in- 
teract freely amongst each other. The unstable particle, 
which was created in interaction process between these 
two theories before taking the limit, becomes so heavy 
that it can not be formed anymore at any energy scale. 



III. FORM FACTORS 

We are now in the position to compute the n-particle 
form factors related to this model, i.e. the matrix el- 
ements of a local operator 0(x) located at the origin 
between a multiparticle in-state of particles (solitons) of 
species [A, created by V fi (9), and the vacuum 

pO^-Mn . . . 9n) = (O(o)Vk {6i)V„ (9 2 ) ...V„ n (9 n )) . 

(5) 

We proceed in the usual fashion by solving the form fac- 
tor consistency equations HQ. For this purpose we ex- 
tract explicitly, according to standard procedure, the sin- 
gularity structure. Since no stable bound states may be 
formed during the scattering of two stable particles the 
only poles present are the ones associated to the kine- 
matic residue equations, that is a first order pole for 
particles of the same type whose rapidities differ by iir. 
Therefore, we parameterize the n-particle form factors as 
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As usual we abbreviate the rapidity difference as 0jj = 
9i — 9j. Aiming towards a universally applicable and 
concise notation, it is convenient to collect the particle 
species . . ji n in form of particular sets 



Wk{k) = {n\n = i} 
m ± (i ± )= \Jm(k) 



(7) 
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The number of elements belonging to the sets 97ti,9Jt± 
is indicated by their arguments k,l±, respectively. We 
understand here that inside the sets 9Jl± the order of the 
individual sets OJli is arbitrary. This simply reflects the 
fact that particles of different species but identical colour 
interact freely. However, DJl is an ordered set since ele- 
ments of 0JI+ and 97t_ do not interact freely and w.l.g. 
we adopt the convention that particles belonging to the 
"+" -colour set come first. The H n are some overall con- 
stants and the Q n are polynomial functions depending on 
the variables Xi = exp#i which are collected in the sets 
X,Xi,X± in a one-to-one fashion with respect to the parti- 
cle species sets DJl,Wl\ ,97t± . The functions (dij ) are 
the so-called minimal form factors which by construc- 
tion contain no singularities in the physical sheet and 
solve Watson's equations ||,|9| for two particles. For the 
SU{N) 2 -HSG model they are found to be 



FZM=N^ (sin- 
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Here hf = 2* exp (wr(l - c;)/4 + c;#/4 - G/4) is a nor- 
malization function with G being the Catalan constant. 
It is convenient also to introduce the function F^ in (9) = 
^ e -ci6/4,pV) (0^y*i . The minimal form factors obey the 
functional identities 
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Substituting the ansatz (^|) into the kinematic residue 
equation ||||[l| , we obtain with the help of ([fl]) a recur- 
sive equation for the overall constants for /Ltj G 9JI+ 



H 
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We introduced here the numbers Zj = ^ 69J1 
which count the elements in the neighbouring sets of 9Jti. 
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The Q-polynomials have to obey the recursive equa- 
tions 



n O\Wl(2+l + ,l-) 
^n+2 



k=0 

x( _ j)2Sl+ r 1 + l CT2s!+ri(Xi)Q 0| OT ( i+ , L ) (x) (13) 

For convenience we defined the sets X xx :— {—x, x} U X, 
x :=ie ai - i+1 X and the integers Q which are or 1 depend- 
ing on whether the sum $ + r, is odd or even, respec- 
tively. $ is related to the factor of local commutativity 
u> = (— 1) = ±1. (7k(xx, . . . , x m ) is the fc-th elemen- 
tary symmetric polynomial. Furthermore, we used the 
sum convention lyXj = U^esm hj^j an d parameterized 
li = 2si +Tj, U = 2Si + ?i in order to distinguish between 
odd and even particle numbers. 

We will now solve the recursive equations ( |l2| ) and 
( |l3| ) systematically. The equations for the constants are 
solved by 



Si(2s i -i i -l-2r i ) 



- H 0\ n ,u 



(14) 

The lowest nonvanishing constants H°^ n ' li are fixed by 
demanding, similarly as in the Si7(3) 2 -case that 
any form factor which involves only one particle type 
should correspond to a form factor of the thermally per- 
turbed Ising model. To achieve this we exploit the ambi- 
guity present in (|l2|), that is the fact that we can multiply 
it by any constant which only depends on the /--quantum 
numbers. 

As the main building blocks for the construction of the 
Q-polynomials serve the determinants of the (t + s) x(t + 
s) -matrix 



o-20-i) + „+(£ + ),l <i<t 

(T 2 {j-i+t) + v- (£-) ,t<i<s + t 

(15) 

for , t~^~ — 0, 1 which were introduced in Q]. The deter- 
minant of A essentially captures the summation in fll3| ) 
due to the fact that it satisfies the recursive equations 



det A 



2+Z,2t+T- ^ + 



\P=0 / 

xdet^^;_(X+,i_) (16) 

as was shown in Analogously to the procedure in 
Q we can build up a simple product from elementary 
symmetric polynomials which takes care of the pre-factor 
in the recursive equation (|l3|). Defining therefore the 
polynomials 

q^-)(x + ,x_)= n^^w 3 ^ 7 ^ (i7) 

p i/fc eOT +/ _ 

X(T 2si+n (Xj) 2 Oj \I,,X,) 2 Ol k (X k )) 2 



it follows immediately with the help of property (|T(]) that 
they obey the recursive equations 



-,9JI(2+Z + ,Z_) /«i ™ \ 
S2+n V x + ' x - I 

Si 

^{si- P )+n+i-u 



f( / + J -)(X + ,X_)a 2Si+ri (X. i ) 
a 2p +^(IijXj) . (18) 



x^x 

p=0 



Comparing now the equations (|13|) and ( pL 8[ ) we obtain 
complete agreement. Notice that the numbers T>i are not 
constrained at all at this point of the construction. How- 
ever, by demanding relativistic invariance, which on the 
other hand means that the overall power in ([j|) has to be 
zero, we obtain the additional constraints 



Vi = 1 + Tj - V l 



and 



= r^Vi - 1) . (19) 



Taking in addition the constraints into account which are 
needed to derive ( |l6|) (see Q), this is most conveniently 
written as 



2 + Ti > n, 



i>i > n 



(20) 



For each fii G 97t+ the equations (|2(]) admit the 10 fea- 
sible solutions found in [pi. However, one should notice 
that the individual solutions for different values of i are 
not all independent of each other. We would like to stress 
that despite the fact that ( p7| ) represents a huge class of 
independent solutions, it does certainly not exhaust all 
of them. Nonetheless, many additional solutions, like the 
energy momentum tensor, may be constructed from (|2(i| ) 
by simple manipulations like the multiplication of some 
CDD-like ambiguity factors or by setting some expres- 
sions to zero on the base of asymptotic considerations 
(see ||) for more details. For many applications we wish 
to carry out in the next section we require the form fac- 
tors for the trace of the energy momentum 0, the first 
non- vanishing terms read 



-2mm 1 sinh(6>/2) 
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i<3 



-2cosh(6»i 2 /2)cosh(6» 3 4/2) 
2 (3+Ecosh(%))n^'( 
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4IL 
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cosh(%/2) 



F, 



J (3+Ecosh(^))n^ 



i<3 



1<J 



4cosh(6»i 2 /2)cosh(6» 3 4/2) 



for lij ^ and Ikj ^ . When considering the RG flow 
in the next section, it will be important to note that from 



lim. 



F; 



(6>) ~ exp(— CT^ i+ i/4) follows 



lim irflw+i-(0) =0 

<Ji i-fl— *O0 



(21) 



Having determined the form factors, we are in princi- 
ple in the position to compute the two-point correlation 
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function between two local operators in the usual way by 
expanding it in terms of n-particle form factors 



oo - 

<0(r)0'(O))=£ E J 



d9\ . . . d6 n ^-rE 



n—1 . ■■fx n _ 



n\(2n) n 



(22) 



We abbreviated the sum of the on-shell energies as E = 
yV, m Ui cosh^. Now we want to evaluate the expres- 
sion ( p2| ) in several different applications in order to com- 
pute various quantities of interest. 



IV. RENORMALIZATION GROUP FLOW 

Renormalization group methods have been developed 
originally JTo) to carry out qualitative analysis of regions 
of quantum field theories which are not accessible by per- 
turbation theory in the coupling constant. In particular 
the /3-function provides an inside into various possible 
asymptotic behaviours and especially it allows to identify 
the fixed points of the theory. We now want to employ 
this method to check our solutions and at the same time 
the physical picture advocated for the HSG-models. 

For this purpose we want to investigate first of all the 
renormalization group flow, in a similar spirit as for the 
SU (3)2/U (l) 2 -case in ||, by evaluating the c-theorem 

fit 



c(r ) 



drr 3 (e(r)6(0)) 



(23) 



In particular for r$ — the function c(ro) coincides with 
Ac = c uv — q r , i.e. the difference between the ultra- 
violet and infrared Virasoro central charges. Comput- 
ing the correlation function for the trace of the energy- 
momentum tensor in ( |23| ) by means of ( p2[ ) and using 
the form factor expressions of the previous section the 
individual n-particle contributions turn out to be 



Ac (2) = (N — 1) - 0.5 

Ac (4) = (N - 2) ■ 0.197 

Ac (6) = (N - 2) ■ 0.002 + (N - 3) • 0.0924 



5> 

fc=2 



N * 0.7914 - 1.1752 



(24) 
(25) 
(26) 

(27) 



Apart from the two particle contribution (pi[), which is 
usually quite trivial and in this situation can even be 
evaluated analytically, we have carried out the multidi- 
mensional integrals in ( p2| ) by means of a Monte Carlo 
method. We use this method up to a precision which is 
higher than the last digit we quote. For convenience we 
report some explicit numbers in table 1. 
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c 


Ac™ 


Ac< 4 ) 


Ac( 6 > 




3 


1.2 


1 


0.197 


0.002 


1.199 


4 


2 


1.5 


0.394 


0.096 


1.990 


5 


2.857 


2 


0.591 


0.191 


2.782 


6 


3.75 


2.5 


0.788 


0.285 


3.573 


7 


4.6 


3 


0.985 


0.380 


4.365 


8 


5.6 


3.5 


1.182 


0.474 


5.156 



Table 1: n-particle contributions to the c-theorem versus 
SU(N) 2 /U(1) N - 1 -WZNW coset model central charge. 

The evaluation of (p4])-(p7|) illustrates that the scries 
( p2| ) converges slower and slower for increasing values of 
N, such that the higher n-particle contributions become 
more and more important to achieve high accuracy. Our 
analysis suggests that it is not the functional dependence 
of the individual form factors which is responsible for 
this behaviour. Instead this effect is simply due to the 
fact that the symmetry factor, that is the sum ^ , 
resulting from permutations of the particle species in- 
creases drastically for larger N. 

Having confirmed the expected ultraviolet central 
charge, we now study the RG-flow by varying tq in (|23|) . 
We expect to find that whenever we reach an energy scale 
at which an unstable particle can be formed, the model 
will flow to a different coset. This means following the 
flow with increasing tq we will encounter a situation in 
which certain Cj^+i are considered to be large and we ob- 
serve the decoupling into two freely interacting systems 
in the way described in (||). For instance for the situ- 
ation 012 > (T23 > CT34 > ... we observe the following 
decoupling along the flow with increasing tq: 

SU(N) 2 /U(1) N - 1 

I 

SU(N - l) 2 /U(l) N - 2 ® SU(2) 2 /U(1) 
I 

SU(N - 2) 2 /U(l) N - 3 ® (SU(2) 2 /U(\)f 

i 



I 

(SU(2) 2 /U(1)) 



N-l 



We can understand this type of behaviour in a semi- 
analytical way. The precise difference between the central 
charges related to (^ is 



c su{i+i) 2 /u(iy®su(N~i) 2 /u(i) N -'- 1 - 

2i(N + 5)(N 



(28) 



c SC/(7V) 2 /;7(l) JV - 1 



(N + 2)(i + 3)(N -i + 2)' 
Noting with fl21|) that we loose at each step all the con- 



tributionsFn to Ac, we may collect the val- 



ues (|24|)-(26), which we have determined numerically and 
find 



lim 

i+i— > 



•) 



Ac Ki+i = 0, • • ■) - 0.2914I M+1 -0.0924I M _i 



(29) 
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Similarly as for the deep UV-region we find a relatively 
good agreement between (^8|) and ((2^) for small values 
of N. The difference for larger values is once again due 
to the convergence behaviour of the series in (|22]). 

For ro — qualitatively a similar kind of behaviour 
was previously observed, for the two particle contribution 
only in the context of the roaming Sinh-Gordon model 
|jl2"|| . Nonetheless, there is a slight difference between 
the two situations. Instead of a decoupling into different 
cosets in these type of models the entire S-matrix takes 
on the value —1, when the resonance parameter goes to 
infinity. The resulting effect, i.e. a depletion of Ac, is 
the same. However, we do not comply with the interpre- 
tation put forward in p2| , namely that such a behaviour 
should constitute a "violation of the c-theorem". The 
observed effect is precisely what one expects from the 
physical point of view and the c-theorem. 

We present our full numerical results in figure 1 , which 
confirm the outlined flow for various values of N. 




To establish more clearly that the plateaus admit in- 
deed an interpretation as fixed points and extract the def- 
inite values of the corres pon ding; Virasoro central charge 
we can also, following determine a (3 type func- 



tion from c(r). The 
Symanzik equation JU 



function should obey the Callan- 



r Tr 9 = PQ,) 



(30) 



The "coupling constant" g :— c uv — c(r) is normalized 
in such a way that it vanishes at the ultraviolet fixed 
point. Whenever we find (3(g) — 0, we can identify c = 
c U v — g as the Virasoro central charge of the corresponding 
conformal field theory. Hence, taking the data obtained 
from (p3|), we compute f3 as a function of g by means of 
d30| ) . Our results for various values of N are depicted in 
figure 2, which allow a definite identification of the fixed 
points corresponding to the coset models expected from 
the decoupling (0). 



-40 -30 

t„ = 2 In r„ 
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Figure 1: RG flow for the Virasoro central charge. 

We observe that the c-function remains constant, at a 
value corresponding to the new coset, in some finite inter- 
val of ro. In particular, we observe the non-equivalence of 
the flows when the relative order of magnitude amongst 
the different resonance parameters is changed. For N = 5 
we confirm (we omit here the ?7(l)-factors and report the 
corresponding central charges as superscripts on the last 
factor) 



Pi 2 > g23l SU(5) 2 7 g23 > trial 

y \ 

SU(A) 2 (g)SU(2)l SU(3) 2 <8> SU{3)f 

\ / 
SU(3) 2 <8> SU{2) 2 <g> SU(2)J 
I 

SU(2) 2 <gi 517(2)2 <8 SU{2) 2 eg) 517(2)3 

The precise difference in the central charges is ex- 
plained with (|2g|), since the contribution 0. 09247;. j_i only 
occurs for i = 1. 



Figure 2: The /3-function. 

For SU(A) 2 we clearly identify from figure 2 the four 
fixed points g — 0,0.3,0.5,2 with high accuracy. The 
five fixed points g = 0,0.357,0.657,0.857,2.857, which 
we expect to find for SU(5) 2 are all slightly shifted due 
to the absence of the higher order contributions. 



V. OPERATOR CONTENT OF S , (7(V) 2 /C/(l) iv ~ :L 

We now want to identify the operator content of our 
theory by carrying out the ultraviolet limit and match- 
ing the conformal dimension of each operator with the 
one in the SU(N) 2 /U(1 ) Ar - 1 -WZNW-coset model. For 
this purpose we have to determine first of all the entire 
operator content of the conformal field theory. 

According to |T^] the conformal dimensions of the 
parafermionic vertex operators are given by 



A(A,A) 



(A-(A + 2p)) (A -A) 
(4 + 2 N) 4 



(31) 
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Here A is a highest dominant weight of level smaller or 
equal to 2 and p = l/2^ a >o a * s * ne Weyl vector, i.e. 
half the sum of all positive roots. The A's are the cor- 
responding lower weights, which may be constructed in 
the usual fashion (see e.g. llq]): Consider a complete 
weight string A + na, . . . , A, . . . , A — ma, that is all the 
weights obtained by successive additions (subtractions) 
of a root a from the weight A, such that A + (n + 1)q 
(A — (m + l)a) is not a weight anymore. It is then a well 
known fact that the difference between the two integers 
to, ft is m — n = X ■ a for simply laced Lie algebras. This 
means starting with the highest weight A, we can work 
our way downwards by deciding after each subtraction 
of a simple root cc, whether the new vector, say x, is a 
weight or not from the criterion rrii = rii + x ' a i > 0- 
With the procedure just outlined we obtain all possible 
weights of the theory. Nonetheless, it may happen that a 
weight corresponds to more than one linear independent 
weight vector, such that the weight space may be more 
than one dimensional. The dimension of each weight vec- 



tor is computed by means of 



£ a >oE^l 2 "A+j a ((A + fa)-Q) 

((A + A + 2p)-(A-A)) 



(32) 



For consistency it is useful to compare the sum of all these 
multiplicities with the dimension of the highest weight 
representation computed directly from the Weyl dimen- 
sionality formula (see e.g. |l5[) 



E 



= dim A = 



n 

Q>0 



((A + p)-a) 
(p-a) 



(33) 



To compute all the conformal dimensions A (A, A) accord- 
ing to (3l|) in general is a formidable task and therefore 
we concentrate on a few distinct ones for generic N and 
only compute the entire content for N = 4. 

Noting that Ai • Xj — K~j X , with K being the Cartan 
matrix, we can obtain relative concrete formulae from 
(|3l|). For instance 



A(A l ,A l ) 



JV-1 i 



8 + 2N 



(34) 



Similarly we may compute A(Ai + Xj, Xi + Xj), etc. in 
terms of components of the inverse Cartan matrix. Even 
more explicit formulae are obtainable when we express 
the simple roots a* and fundamental weights A^ of SU (N) 
in terms of a concrete basis. For instance we may choose 
an orthonormal basis {ei} in M. N (see e.g. Jl6)), i.e. Si ■ 



i . N 

•et+i, Xj = y^gj - — y^gj, % 



.N-l. 



Noting further that the set of positive roots is given by 
|gj — £j : : 1 < i < j ' < N}, we can evaluate (§l]), @) and 
(|33|) explicitly. This way we obtain for instance 



A(AiA) = y and A(2\, 2A 4 ) = 0. (35) 

o + 4iV 

Of special physical interest is the dimension of the per- 
turbing operator. As was already argued in it cor- 
responds to A(ip,0), with t/j being the highest root, and 
moreover it is unique. Noting that for SU(N) we have 
ip = Ai + Ajv— i, we confirm once more 



A(^ = Ax+Ajv-x.O) = 



N 



(36) 



A^ + 2 

Other dimensions may be computed similarly. 

A. The SU{4)/U{lf example 

For SU(4)/U(1) 3 we present the result of the compu- 
tation of the entire operator content in table 2. In case 
the multiplicity of a weight vector is bigger than one, we 
indicate this by a superscript on the conformal dimen- 
sion. 



A\A 


Xi 


A 2 


2Ai 


2A 2 


Ai + A 2 


Ai + A 3 


dim A 


4 


6 


10 


20 


20 


15 


A 


1/8 


1/6 








1/8 


1/6 


A — Qi 


1/8 




1/2 




1/8 


1/6 


A — Q2 




1/6 




1/2 


1/8 




A — a 3 












1/6 


A — Oil — OL2 


1/8 


1/6 


1/2 


1/2 


5/8" 


1/6 


A — «2 — «3 




1/6 




1/2 


1/8 


1/6 


A — ol\ — as 












1/6 


A - 2ai 















A - 2a 2 















A - 2ai - 2a 2 
















A - 2a 2 - 2a 3 















A — 2qi — a?2 






1/2 




1/8 




A — ai — 2a2 








1/2 


1/8 




A — 2q2 — us 








1/2 






A — a?i — a.2 — Q3 


1/8 


1/6 


1/2 


1/2 


5/8" 


2/3 a 


A — 2cn — a.2 — ot3 






1/2 




1/8 


1/6 


A — ai — 2a2 — «3 




1/6 




i 2 


5/8" 


1/6 


A — Oil — OL2 — 2o>3 












1/6 


A — 2ol\ — 2ct2 — Q!3 






1/2 


1/2 


5/8" 


1/6 


A — ai — 2a2 — 2a3 








1/2 


1/8 


1/6 


A - 2ai - 2a 2 










1/8 




A - 2ai - 2a 2 - 2a 3 












1/8 


1/6 


A — qi — 3a2 — 2a3 








1/2 






A — on — 3a2 — ct'j, 








1/2 






A — 2ol\ — 3Q2 — Ct3 








1/2 


1/8 




A — 2ai — 3o?2 — 2«3 








1/2 


1/8 




A - 2ai - 4a 2 — 2a 3 
















Table 2: Conformal dimensions for O a( - a ' x) in the 
S77(4)2/(7(l) 3 -WZNW coset model. 

The remaining dominant weights of level smaller or 
equal to 2, namely A = A3, 2 A3, A2 + A 3 , including their 
multiplicities may be obtained from table 2 simply by the 



G 



exchange 1 <-> 3, which corresponds to the Z2-symmetry 
of the S[/(4)-Dynkin diagram. 

Summing up all the fields corresponding to different 
lower weights, i.e. not counting the multiplicities, we 
have the following operator content 

2/3 ,0\u x O ,8 x C 5/8 ,18 x C 1/6 ,24 x C9 1/2 ,32 x C 1/8 , 
that is 98 fields. 

VI. OPERATOR CONTENT OF HSG 

We will now turn to the massive model and evaluate 
the flow of the conformal dimension |3) 



A°(r ) 



2(0(0)) 



drr (0(r)O(O)> 



(37) 



Here O is a local operator which in the conformal limit 
corresponds to a primary field in the sense of p7[ . In 
particular for ro = 0, the expression ( |37|) constitutes the 
delta sum rule [p!s[ , which expresses the difference be- 
tween the ultraviolet and infrared conformal dimension 
of the operator O. 

We start by investigating the operator which in the 
case when all particles in (||) are of the same type cor- 
responds to the disorder operator /i in the Ising model. 
Using the fact that we should always be able to reduce 
to that situation, we consider the solution corresponding 
to n — fi — v L = ft = for all i. Then the A-sum rule 
( |37j ) yields for the individual n-particle contributions 



A m(2) = (tv - 1) • 0.0625 
A m(4) = (2 - JV) ■ 0.0263 
A^ 6) = (N - 2) ■ 0.0017 



k=2 



0.0266 



(3 

N * 0.0206 . 



(38) 
(39) 

N) ■ 0.0113 (40) 
(41) 



We assume that this solution has the conformal dimen- 
sion A(Ai, Ai) in the ultraviolet limit. For comparison we 
report a few explicit numbers in table 3. 



N 


A(A x ,Ai) 


A^( 2 ) 


A^( 4 ) 


AM 6 ) 


Efe=2 M 


3 


0.1 


0.125 


-0.0263 


0.0017 


0.1004 


4 


0.125 


0.1875 


-0.0526 


-0.0079 


0.1270 


5 


0.143 


0.25 


-0.0789 


-0.0175 


0.1536 


6 


0.156 


0.3125 


-0.1052 


-0.0271 


0.1802 


7 


0.16 


0.375 


-0.1315 


-0.0367 


0.2068 


8 


0.175 


0.4375 


-0.1578 


-0.0463 


0.2334 



Table 3: n-particle contributions to the A-theorem versus 
conformal dimensions in the SC/(W) 2 /(7(1) JV ~ 1 -WZNW coset 
model. 

As we already observed for the c-theorem, the series 
converges slower for larger values of N. The reason for 



this behaviour is the same, namely the increasing sym- 
metry factor. Note also that the next contribution is 
negative. 

Following now the RG-flow for the conformal dimen- 
sion (j37j) by varying ro, we assume that the A(Ai, Air- 
field flows to the A(Ai, Ai)-field in the corresponding new 
cosets. Similar as for the Virasoro central charge we may 
compare the exact expression 



A(Ai,Ai) 



SU(i+l) 2 /U(l) i ®SU(N-i) 2 /U(iy 



(42) 



A n w i{N + S)(N-i-l) 

X ^su(N h /u(i)^ + m + m + 3){N _ l + 2 y 



with the numerical results. The contributions (pq)-(EQ) 
yield 

A*K i+1 ,...) = 

0, . . .) + 0.03597, . 



lim 

i+l— * 



(43) 



A M (er 



0.01131; 



Once again we find good agreement between the two com- 
putations for small values of N. Our complete numerical 
results are presented in figure 3, which confirm the out- 
lined flow for various values of N. 
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Figure 3: RG flow for the conformal dimension of /x. 

Notice by comparing the figures 3 and 1, that, as we 
expect, the transition from one value for A to the one in 
the decoupled system occurs at the same energy scale to 
at which the value of the Virasoro central charge flows to 
the new one. 

In analogy to (|3(]) we may now define a function "/3"' 
and demand that it obeys the Callan-Symanzik equation 

r|-s' = W) ■ (44) 
ar 

The "coupling constant" related to f3' is normalized 
in such a way that it vanishes at the ultraviolet fixed 
point, i.e. g' :— A(r) — A uv , such that whenever we 
find P'(g') — 0, we can identify A = g' — A uv as the 
conformal dimension of the operator under consideration 
of the corresponding conformal field theory. From our 
analysis of ( p7|) w e may determine j3' as a function of g' 
by means of (E4I). Our results are presented in figure 4. 



7 




-0.05 0.00 
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Figure 4: The /^'-function . 

Once again, for SU (4)2 the accuracy is very high and 
we clearly read off from figure 4 the expected fixed points 
g' = -0.125,0,0.0375,0.0625. The ST/(5) 2 -fixed points 
g' = -0.1429,0,0.0446,0.0821,0.1071, are once again 
slightly shifted. 
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Figure 5: Rescaled correlation function G(R) = (B(i?)0(O)) as 
a function of R = rm. 

Unfortunately, whenever the correlation function be- 
tween O and © is vanishing, or when we consider an 



operator which does not flow to a primary field, we can 
not employ the delta sum rule (|37|). Alternatively, we 
may exploit the well known relation 



Hm (O(r)O(O)) 

r— vu 



-4 A" 



(45) 



near the critical point in order to determine the confor- 
mal dimension. To achieve consistency with the proposed 
physical picture we want to identify in particular the con- 
formal dimension of the perturbing operator. Recalling 
that the trace of the energy momentum tensor is propor- 
tional to the perturbing field we analyse (0(i?)0(O)) for 
this purpose. 

According to (^5|), we deduce from figure 5 A = 
2/3,5/7 for N — 4, 5, respectively, which coincides with 
the expected values. 



VII. CONCLUSIONS 

One of the main deductions from our analysis is that 
the scattering matrix proposed in || may certainly be 
associated to the perturbed gauged WZNW-coset. This 
is based on the fact that we reproduce all the predicted 
features of this picture, namely the expected ultraviolet 
Virasoro central charge, various conformal dimensions of 
local operators and the characteristics of the unstable 
particle spectrum. 

Our construction of general solutions to the form factor 
consistency equations certainly constitutes a further im- 
portant step towards a generic group theoretical under- 
standing of the rt-particle form factor expressions. The 
next natural step is to extend the investigation towards 
higher level algebras p0[ . 

Concerning the computation of correlation functions, 
our results also indicate that the "folkloristic belief" of 
the fast convergence of the series expansion of ( ^2| ) has 
to be challenged. In fact, for large values of N, this is 
not true anymore. It would be highly desirable to have 
more concrete quantitative criteria at hand. 

Despite the fact of having identified some part of the 
operator content, it remains a challenge to perform a def- 
inite one-to-one identification between the solutions to 
the form factor consistency equations and the local oper- 
ators. It is clear that we require new additional technical 
tools to do this, since the A-sum rule ( |37j ) may not be 
applied in all situations and ( f45| ) does not allow a clear 
cut deduction of A. 

In comparison with other methods to achieve the same 
goal, we should note that in principle we could obtain, 
apart from conformal dimensions different from the one 
of the perturbing operator, the same qualitative picture 
from a TBA-analysis [^9|. However, in the latter ap- 
proach the number of coupled non-linear integral equa- 
tions to be solved increases with N, which means the 
system becomes extremely complex and cumbersome to 
solve even numerically. Computing the scaling function 
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with the help of form factors only adds more terms to 
each n-particle contribution, but is technically not more 
involved. The price we pay in this setting is, however, 



the slow convergence of 



We conjecture that the "cutting rule" which describes 
the renormalization group flow also holds for other 
groups. This is supported by the general structure of 
the HSG-scattering matrix. 
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